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I - Introduction 



Several strongly correlated fermionic systems such as liquid He, heavy- fermion com- 
pounds, high-T c superconductors, and Kondo systems are the subject of intense theoretical 
as well as experimental interest. The Hubbard model, originally introduced to describe 
correlation effects in narrow <i-band materials, has been put forward as a possible key to 
the understanding of high-T c oxide superconductivity [1] . 

Describing electrons on a lattice with one orbital per site, the Hubbard Hamiltonian : 

H = J2 Uj c\ a c 3a + ni*ni(-a) (1) 

ij,a i 

includes a kinetic hopping term tij between different sites and an on-site Coulomb repul- 
sion U between electrons of different spin. The operator Ci„ annihilates an electron of 
spin a at site i and the occupation number operator is hi a = c\ a di a . 

Of particular interest is the strong-coupling regime where an effective Hamiltonian of 
the Hubbard model with nearest-neighbor hopping is the t-J model Hamiltonian [2] : 

ff = <E (4a> + h.c.) + j E • ^ - fa ( 2 ) 

<ij>,a <ij> 

( with J = At 2 /U and S- = ±^ cJ a a a/3 c l/3 ) 

a/3 

The on-site Coulomb repulsion there is very large as compared with the electron-hopping 
energy, and therefore when there is less than half filling the system will avoid configuration 
with doubly occupied sites. One has thus the constraint : 

E 6 Ua< 1 (3) 
u 

Apart from numerical approaches, a popular analytical approach to the t-J model 
is the slave-particle theory [3] where for instance in the slave-fermion representation the 
electron operator c\ a is written as c\ a = b\ a fi , with / the slave fermion and b a a 
boson. Instead of the constraint of Eq. (3) , which is difficult to handle, one considers the 
a priori more convenient slave-particle constraint avoiding double occupancy at site i : 

/M + E^- = 1 ( 4 ) 

a 

With the boson operator b\ a keeping track of the spin and the fermion operator fj 
keeping track of the charge, this formalism is well adapted to study the problem of the 
decoupling of spin and charge degrees of freedom in the large- U- limit Hubbard model. This 
decoupling, characteristic of Luttinger liquids, appears to occur in ID [4], the situation 
being still confused in 2D. At the mean field level, spinons and holons, the elementary 
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spin and charge excitations, may be separated in this formalism but are strongly coupled 
beyond the mean-field approximation. 

Such slave-particle approaches are usually studied in a functional integral (over co- 
herent states of Fermi and Bose fields) representation of the partition function, the slave- 
particle constraints being enforced by functional integration over Lagrange multipliers. 

However one may wonder what happens within a direct operator quantum approach 
of such slave-particle theories. Indeed, sum rules, coming from operator commutation 
relations, for spectral functions of the boson and fermion were used recently in Ref. [5] in 
a study of the electron-momentum distribution function in the t-J model in the framework 
of the slave-particle approach and within the decoupling scheme for the electron Green's 
function. It was claimed there that the sum rule for the physical electron was not obeyed 
within this framework and correspondingly that the electron Fermi surface (EFS) was not 
explained. However it happens that the operator commutation relations used in Ref. [5] for 
the slave fermion and boson were the usual naive ones, i.e. the same as if the slave- particle 
constraint was not present. One then can question the use of such operator commutation 
relations in a slave-particle approach. 

In this paper, we study the direct operator approach of such slave-particle theories. In 
Sec. II, we present at the classical level the consistent Hamiltonian formulation of models 
having a slave-particle constraint for their fields. We show in Sec. Ill at the quantum level, 
for the slave- fermion and the slave-boson representations of the t-J model, the modifications 
in the sum rules for the slave particles coming from the fact that the correct canonical 
relations compatible with the constraints are not the naive ones. We present in Sec. IV 
a direct explicit operator quantization of the slave-particle approaches of the t-J model 
which confirms the results obtained in Sec. Ill . In Sec. V, we extend our analysis to a 
slave-boson representation which has been introduced [6] for the finite- U Hubbard model. 
Sec. VI summarizes our conclusions, and an Appendix presents calculations omitted for 
clarity in Sec. II. 

II - Hamiltonian formulation for slave models 

Let us consider the general classical Lagrangian (written in real time) for n bosons 
b a and m fermions f a (Grassmann variables) on a lattice : 

L = i E b l9Aa + i £ fldtU + X{b\ 6, /t, /) + A, ( E b l b ™ + E fUur - 1 ) 

i,a i,a i a a 

(5) 

bosons and fermions being submitted to the slave- particle constraint at each site % : 

*i = + -1 = (6) 

a a 

One shows in the Appendix that, after an d la Dime treatment [7], one gets an 
Hamiltonian formalism with the non zero brackets : 

\bia, b )p] = ~ i <%<W {Via, /J T } = - i SijSaT {A l7 II j} = Sij (7) 
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but with the first class [7-9] constraints 



$i«0 and Ui « (8) 

where, following Dirac, first class means here {$^,11^} tv and where the symbol ~ 
(weakly zero) means that one has to set the constraints only after computing all the 
brackets. We use graded Poisson-Dirac brackets [9] such that for instance : 

{B, A} = - (-l) ab {A, B} {A, BC} = {A, B}C+ (-l) ab B {A, C} (9) 

where a = if A is a bosonic quantity and a = 1 if A is a fermionic quantity. 

The role of the first class constraints is to generate infinitesimal contact transforma- 
tions (that we shall call as usual [8] gauge transformations) in the Hamiltonian formalism 
that do not affect the physical state of the system. We have here for the non-zero brackets 
of the fundamental variables with the first class constraints : 

V>ia, *j} = ~ * haSij [b\ a , = + % bjjij (10) 

{f ia , = - i f ta 5 tJ *,-} = + i fl6 tj (11) 

{A i ,U j } = S ij (12) 

which means that the gauge transformations are b — >■ e~ t0 b , — > e l9 tf , / — > e~ ld f , — > 
e l9 p ,A^A + a,n->n . 

The standard strategy [8] is then to fix the gauge by choosing explicit forms for each 
gauge and imposing them as constraints not following from the Lagrangian. The choice of 
gauges should be made in such a way that the constraints <E>j and Ilj will cease to be 
first class. It happens that a convenient choice here is the linear one : 

d>; = ( G ia b ia + G\ a b\ a ) + ( Hurfi* ~ H\jl ) + K t « (13) 

a a 

A, « (14) 

where the G 's, the K 's and the (Grassmannian) H 's are parameters. 

We have then for the non-zero brackets of the fundamental variables with these gauge 
fixing constraints : 

{b ia , *;} = -i G\j tJ {b\ a , $;•} = + i Gi a Sij (is) 

*J} = - i fT^tfy {/t + i ff^y (16) 

{U i ,A j } = -S ij (17) 

and we obtain the bracket of the slave-particle constraint with its gauge fixing constraint 
as : 



{*i,*;} = i 



^ ( G ia b ia - G\ a b\ a ) + H ia f ia + H\j\ a ) 



5ij = i Di Sij (18) 
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which is not weakly zero. 

Let us note that the Hamiltonian (see the Appendix) is now : 

h 3 = - x{b\ b, /t, (xi - Ai) ( b l b ™ + E flf* - 1 ) + E u ^ ( 19 ) 



where the first class constraints coefficients x and w (which are in H3 independent of the 
fields) are determined by the requirement that the time derivative of the gauge fixing 
constraints is weakly zero : 

r) X f)X f)X f)X 

^3} = i E ( G ^7-r - G l^-) E - H lar) - i & - A *) D * « 

(20) 

{Ai,H 3 } = Wi ^0 (21) 

Defining ^1 = *i > ^2 = , ^3 = ^ , y? 4 = A; , the matrix C a6 = {<^ a ,v3 6 } 
is non singular and all the constraints are now second class [7-9]. Systematic use of the 
standard Dirac bracket [7] of two quantities A and B : 

{A, B}^ = {A, B}-J2 {A Va} (C-% Wb, B} (22) 

a, b 

then allows one to set all these second class constraints strongly to zero because the Dirac 
bracket of anything with a second class constraint vanishes. 

We thus obtain the correct classical non zero canonical relations for slave-particle 
models, compatible with the constraints : 

i { b L b) p }^ = [ ( G ifl bl - G ia bl )/D i ] S tJ (24) 
i {Ka, b jP }«= [ ( G\ p b ia - G\ a b ip )/D i } (25) 

1 { fi ^ = [ 6aT + ( Htrfia + H * f * ),Di ] Sij (26) 

1 {/i ' 4}, = " t ( + H ia fl )/D i ] tf y (27) 

i {/*r> /ir}* = - [ ( #Jr/ia + )/A ] (28) 

* 6 J«}, = " [ ( G iafia ~ )l D i 1 % (29) 

4 { 6 ;«}„ = " t ( - #*Aa )/A ] 6^ (30) 

< { = [ ( - H ia b\ a )/D i ] ^ (31) 

6 i«h = [ ( G\ a f ia - Hlb ia )/D i } Sij (32) 
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It is then clear, as we shall explicitly show below, that a correct operator quantization of 
such slave-particle models must be based on these new canonical relations, and not on the 
naive ones of Eq. (7) which are not compatible with the constraints. 

Let us note that we have {Ilj, A^}^ = {IL^IIj}^ = {Aj, A^}^ = indeed compatible 
with the constraints Hi = Aj = . On the other hand the Hamiltonian is finally : 

H± = -X{b\bj\f) (33) 

and we have consistently for example {bi a , H 3 } = {&i a , -££4}* since : 

dX 

{b ia , H 3 } = i—^ i(xi- A t )b ice (34) 



{ha, H,}, = l-^j- - I — 



v ax , ax v ax t ax 



/ ItT 



35) 



which are indeed equal using Eq. (20). 



One may wonder at this stage about what happens for the classical version of the 
Sl(l 1 2) superalgebra obeyed in the t-J model [10] by the Hubbard [11] operators Xf b : 

[Xf, Xf] ± = (Xf % c ± Xfd ad ) S iS with X™ + Xr = 1 (36) 

a 

which opens the way for a supersymmetric t-J model [10] and which is verified in the 
literature within a slave representation using the naive canonical relations. In fact, taking 
for instance in the t-J model the slave-fermion representation c\ a = b\ a fi of the electron 
field c\ a which corresponds to X"° , the 57(1 12) superalgebra is obeyed using either 
the naive initial brackets of Eq. (7) or the Dirac brackets of Eqs. (23-32). For example 
one has : 

i {4a> Cjf)} = i {4a> Cjfi}^ = [ 4^ + fjfi 5 a/3 ] 8ij (37) 

where b\ a bip corresponds to Xf^ and fjfi to X®°. The reason for this property 
is that the Dirac bracket of any two gauge invariant quantities is the same as their ini- 
tial bracket. As all the generators of the £7(1 12) superalgebra are gauge invariant, it is 
licit to use the initial brackets of the &'s and fs to compute the Dirac brackets of these 
generators. One must however realize that this property is valid only for gauge invariant 
quantities. At the quantum level, a proper quantization should inherit the same prop- 
erty for (anti) commutators of gauge invariant quantities, but not for products of these 
quantities. Let us also note that, using the slave-particle constraint, one has from Eq.(37): 

{ c L c j a } = E* { c Ls«h = [ 1 + f * fi ] Sij (38) 

a a 

a relation that we shall recover below in a sum rule at the quantum level. 
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Ill - Modifications of the naive sum rules for the slave t-J models 

As mentioned in the Introduction, it has been stressed in the literature [5] that a study 
of the electron-momentum distribution function in the t-J model in the framework of the 
slave-particle approach and within the decoupling scheme would give rise to a violation of 
the sum rule of electron number. However, to obtain this result, sum rules for the slave 
particles using quantization of the naive relations of Eq. (7) were used. On the contrary 
we shall show in this Section that starting from the correct canonical relations compatible 
with the constraints produces modifications in these sum rules for the slave particles in 
such a way that the sum rule of electron number is indeed obeyed. Though we shall present 
in the next Section a direct explicit quantization which confirms this result, we think that 
it is important for the clarity of the exposition to see here the simple structure of these 
modifications in the sum rules. 

a) The slave-fermion representation 

Let us consider in the t-J model the slave-fermion representation c\ a = b\ a fi of the 
electron operator c\ a , with the slave-particle constraint avoiding double occupancy at 
site % ' 

/7/i + = 1 (39) 

a 

As in Ref. [5], we assume that there is no Bose condensation of spinons i.e. the temperature 
of the system is + . The hole-doping concentration 5 is given by 5 = l^f} fij ■ 

For discussing the electron-momentum distribution function, the Matsubara electron 
Green's function in imaginary time [12,13] : 

Ea(r,r) = - (^(^(^^(r) 6^(0)^(0))) (r = i - j) , (40) 

where (...) means the thermodynamical average, was considered in Ref. [5] within the 
decoupling approximation written on Fourier transform : 

E a (k,u) n ) = -k^2j}^2F(q,Um) B a (q + k,Um + u n ) (41) 
q u m 

F is the Green's function for the slave- fermion / and B a is the Green's function for the 
boson b a . Introducing the Lehmann's spectral representations : 




one easily obtains the electron spectral function as : 

A ea (k,uj) = / ^A f {q,uj')A ba {q + k,u + uj'){n F {u')+n B {u + u')} (43) 
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where np and n# are respectively the Fermi and Bose distribution functions. Using the 
expressions for the numbers of slave- fermions, bosons and electrons in state q : 



w) n ha (q) = 



/+oo 
gn F (u;) A f (q, 
-oo 

/+oo 
^n F (oj)A ea (q,oj) , 
-oo 



+00 



^n B (u)A ba (q,u) 



(44) 



the identity : np (w) [ np (oj 1 ) + = ub(uj + uj') [1 — rip (u/) ] , and the sum 

rules : 

one obtains the sum rule for the electron spectral function : 

»+oo 



2tt 



A ea (k,oj) 



q,ot 



fqi fq 



-I + 



(46) 



and the expression for the number of electrons in state k : 



q,a 



( k ) = Jr ^ n ba{q + k) I \iqj\\ ) - j<^ri ba {q + k)n f {q) 



(47) 



q,a 



Let us make the remark, not mentioned in Ref. [5] , that the last equation can in fact 
be directly derived from the very definition of the decoupling approximation : 



E a (r, r) = B a (r, r) F(-r,-r) 



(48) 



where we used 



T T (b ia (r)b} a (0))) = -B a (r,r) 



2V(/ 4 T (r)/ i (0)))=F(-r,-r) (49) 



Then, with rj -> 0+ [12] , one has : 

n e = I ^ £ a (0, -,?) = 1 ^ 5«(0, F(0, i|) = | E (E ( ) (M) ( 50 ) 



using B a (0, —rj) = — / b ia bi a / an d F(0, rj) = — ( fif}) ■ This gives in Fourier transform: 



ra e (fc) 



fqfq 



(51) 



expression which, using the anticommutator of the fermion, can immediately be written 
in the form of Eq. (47) . 

The problem now is to evaluate the ther mo dynamical average of the commutator 
and anticommutator entering in these expressions. Using the quantization from the Dirac 
brackets : 



HAS}* 



A,B 



(52) 



one obtains from our results Eqs. (23) and (26) of the preceding Section the form of the 
commutator and anticommutator at equal times : 



bia, bij a 



( 1 - Qia ) S. 



fii fj 



-I + 



( 1 + Gi ) S. 



(53) 



where we have introduced the operators defined through the quantizations : 

(Giabia- G\ a b\ a )/D i =e ia -> 6 ia (Hifi+ H\f} )/Di = 0, -> 0,(54) 

Let us note that the definition of Di from Eq. (18) gives by quantization the relation : 

J]0 ia + 0, = 1 (55) 

a 

With n = E p e ipn S p and na = £ p E pa , it follows that : 



= 1 + (So 



= 2 



(56) 
(57) 



Since with the hole-doping concentration 5 one has : 

jj^2 n f( ( l) = S jj^2n ba (q) = 1-6 



(58) 



q,a 



one obtains the following results 



2tt 



A eQ (/c,cu) = 1 + (H 



(59) 
(60) 
(61) 



One then sees clearly the modifications of the sum rules and of the results of Ref. [5] 
coming from the presence of the 's in the correct canonical relations. The main point 



n e (k) = (1-8) ( 1+ (s ) ) - ^^n 6a ((/ + /c) n f (q) 

q,a 

j r J2Mk) = (l-S)(l+ (h ) -5) 
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is that our results show that one indeed obtains the expected result that, if 5 holes are 
introduced into the half-filled system, the total electron number per site would be (1 — 5) , 
instead of the (1 — 5) 2 found in Ref. [5] . 

In fact, from the quantization of Eq. (38) at the end of the preceding Section, the 
sum rule for the electron spectral function must be : 



E [ + °° £ A ^ = ( E fc" °U + ) = 1 + * E = 1 + 6 

a J -°° \ a I k 



(62) 



being unchanged with respect to Ref. [5]. Comparing with Eq. (59), this gives : 

Ho) = S (63) 



From our results, the expected sum rule of the electron number is indeed recovered with 
^ S ) = 5 , since jj J2k n e(k) = (1 — 5) ( 1 + ^ So) — 5 ) = (1 — 5) . However this does 
not guarantee the existence of an EFS within the decoupling scheme in the slave-fermion 
approach of the t-J model, and the arguments of Ref. [5] against an EFS still apply. 

The electron operator anticommutator furthermore shows that one has not only the 
average equality ( ©i) = ( So) =5 but also the operator equality 0j = fjfi : on 



one hand quantization of Eq. (38) gives ^ a Cj a ,c| a = 1 + fjfi , while on the other 

hand one gets by expressing each term of the anticommutator with the boson and fermion 
operators 



Cice 7 c\ a 



= 1 + flfi + E ~ fitih* + //(©i - fjfyfi (64) 

4- ^— 'a 



effectively leading to a consistent quantization expressed by ©^ = fjfc . 
b) The slave-boson representation 

Our analysis of the t-J model in the slave-boson representation c ia = y ia hi of the 
electron operator c\ a , with the slave-particle constraint avoiding double occupancy at 
site i : 

b\k + J2fU^ = 1 ( 65 ) 

follows along the same lines. As in Ref. [5] , we assume that there is no Bose condensation 
of holons. The hole-doping concentration 5 is given by 5 = (h\hi 
The Matsubara electron Green's function in imaginary time : 



Ea(r, r) = - ( T T ( b\{r) f ia (r) &(<)) 6,(0) ) ) (66) 
within the decoupling approximation written on Fourier transform reads : 

E a (k,u n ) = - j r ^2^^2B(q,oj m -uj n ) F a (q + k,uj m ) (67) 



and the electron spectral function is : 

-+oo 



/-too 
^A b (q,oj'-oj)A f<J (q + k,oj') [n F (oj') +n B (u)'-u)] (68) 

Using the identity np{ui) [rip (w') + n B — oS) ] = rip (u/) [ 1 + n B {oo' and the sum 

rules analogous to Eq. (45) , one obtains the sum rule for the electron spectral function : 



q \ a 



t 

(q+k)a 



-I + 



+ jjJ2 n f^ +k ^ ( KK] ) ( 69 ) 



and the expression for the number of electrons in state k 



q,a 



n e( k ) = jf^2n fa (q + k) ( b q ,b\) + jj^n fa {q + k) n b (q) (70) 



q,a 



Again, the last equation can in fact be directly derived from the very definition of the 
decoupling approximation : 



E a {r,r) = - F a {r,r) B{-r, -r) 



(71) 



where we used : 



T T ( L(r) fl (0) )) = -F a (r,r) ( T T ( b\(r) 6,(0) ) ) = - B(-r, -r) (72) 



Then, with 77 — > + , one has : 



1,0- 



i, a 



using F a (0,-77) = + \fLfi<?/ an d -8(0,77) = - ^Mjy • This gives in Fourier 
transform : 



q,CT 



5,0- 



expression which, using the commutator of the boson, can immediately be written in the 
form of Eq. (70) . 

Now, with the quantizations : 
( GA - G\b\ )/D i -> 0, ( H ia f ia + Hjjl )/Di - e,. 



^e l(J + e, = i 

(75) 
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one obtains from our results of the preceding Section the relations at equal times 



= ( 1 - Qi ) 6., 



'J 



= ( 1 + &ia ) 



(76) 



and with the same notations as above 



^ [/(g+fc)<r> f(q+k)a 



= 2 + ( ^ S 0a 



= 3 



(77) 
(78) 



Inserting these results in Eqs. (69) and (70), one finally gets : 



+oo 



(79) 
(80) 
(81) 



Again, one sees clearly the modifications of the sum rules and of the results of Ref. 
[5] coming from the presence of the © 's in the correct canonical relations and that with : 



le(fc) = (1 - <5) ( 1 - (S )) + ^J2 U f^ + k ^ 

q,a 

l^n e (fc) = (l-5)(l- (s ) +5) 



"0 



(82) 



one has the correct sum rules for the electron spectral function 



/+oo 
-oo 



A ea (q, u) = 1 + 5 



(83) 



and for the electron number : 



~k 77-e (^) = 0- ~ ^) t i ns tead of the (1 — 5 2 ) found in Ref. [5] 



(84) 



However this does not guarantee the existence of an EFS within the decoupling scheme in 
the slave-boson approach of the t-J model, and the arguments of Ref. [5] against an EFS 
still apply. 

Furthermore the same arguments as above lead, using the electron operator anticom- 
mutator, to the operator equality 0j = b\bi . 
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IV - Explicit quantization of the slave-particle approaches of the t-J model 

In the preceding Section, it was shown that the expected sum rules for the electron 
spectral function and for the electron number are indeed found in the slave-particle ap- 
proaches of the t-J model due to the fact that ^ ©i^ = ( = $ > the operator 0^ 
being the new term which is present in the canonical relation of the slave particle when 
one quantizes the correct (Dirac) brackets compatible with the constraints. 

We furthermore proved using the electron operator anticommutator that 0^ was 
the slave particle number operator. We shall show in this Section that a direct explicit 
operator quantization of the slave-particle approaches of the t-J model effectively confirms 
this result without invoking the electron operator. 

a) The slave-fermion representation 

Let us first take the slave-fermion case, where we obtain from the results of Section 

II: 

ft, /]] = ( 1 + ©i ) Sij [ka, b} a ] _ = ( 1 - ®ia ) kj (85) 

where (cf. Eqs. (53) and (54)) the 's correspond to the quantizations : 

(Hifi+ Hjfj)/Di - 0, ( G ia b ia - G\ a b\ a )/Di - ia (86) 

However, since A = J2 a ( G iQ b iOL - G] a b\ a ) + ( H t fi + H}fJ ) , the explicit expressions 
of the operators in terms of the operators of the fermion and bosons are a priori not 
obvious. We shall now obtain these expressions through the structure of the Fock space. 

From the slave-particle constraint : 

/7/i + = 1 (87) 

a 

the operators fj fi and b\ a bi a must be particle-number operators satisfying : 

UlhY = flfi and {b\ a b ia f = b\ a b ia (88) 

Now it follows from Eq. (85) that : 

(flfi? = flfi ~ fKflfi - Qi)fi (blka) 2 = blka + bl(blka ~ Qia)ka (89) 

and a consistent quantization is thus : 

0* = flfi O lce = b\ a b ia (90) 

which from Eq. (85) leads to the relations : 

fifl = 1 b ia bl = 1 (91) 
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We will show below that the dimension of the Fock space at site i is in fact infinite (but 
without contradiction with the non double electron-occupancy at site i ). Thus Eqs. (91) 
are consistent with the fact that the particle-number operators are not unity. Let us remark 
that we have consistently, in accordance with Eqs. (55) and (87), : 



(92) 



On the other hand, the fundamental Dirac brackets compatible with the constraints 
given by Eqs. (23-32) lead to the following expressions : 



fyL^a}^ = fi 

{bL ^W}, = - b h 5 ^ + b l e ^ 



(93) 
(94) 
(95) 



Quantizing these expressions with the same ordering of operators gives via Eqs. (90) and 
(91): 



1fib\a 

it b- # 



bL (fifi 
t /it 



Si) = and thus ft b\ a = (96) 

Q ia )=0 and thus kJJ = (97) 

&l/?Mla = h \ a {blpbip ~ + fyp <W and thus b i(3 b\ a = 5 af3 (98) 

Using all our results, one easily verifies that the operators fj fi and b\Jbi a are indeed 
particle- number operators : for example fj fi acting on the fermion state fj |0) at site 
% has eigenvalue 1 and acting on the two boson states b\ a |0) at site i has eigenvalue . 

One can check the consistency of our quantization : for example quantizing the ex- 
pression 

* { b L Ml/j}, = - b\ p <W + e b\ a 

with the same ordering of operators gives : 



(99) 



bLMM 



icti 







= - bL <W + bLbipHa = - bL 5 aP + bL 5 aP = 



(100) 
(101) 



Through an explicit operator quantization of the slave-fermion approach of the t-J 
model we have thus obtained, at site i , : 



fifi = 1 b ia b\p = 5 af3 

or in other words : 



/* bl = 



hJl = 



fit fj 



( 1 + fifi ) Sa 



fii b\ a — b\ a fi S. 



= ( <W - blpb ia ) 5ij 

fi bia $ij 



(102) 

(103) 
(104) 
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Looking in the same way as above at the quantization of the classical canonical relations 
involving only creators or annihilators, we found either no informations or identities. Such 
relations might thus be only identities at the quantum level. The results of Eq. (102) will 
however be sufficient for the purposes of this paper. 

In fact, as a consequence of this study we directly recover the expressions : 

ft. 



Qi = flfi (ei) = (Eo) = 6 (105) 

which through the analysis of Section III leads to the conclusion that the expected sum 
rules for the electron spectral function and for the electron number are indeed found in 
the slave-fermion approach of the t-J model within the decoupling approximation. Using 
Eq. (50) and our result fifj = 1 , one can also see directly that : 

= vji*) ) ( /*/?) = a - s) ( 106 ) 



It is also important, with the expression of the electron operator 



c, 



t 



bl fi (107) 



ia ia J 1 



and using our results for the quantization, to found the exact electron number operator : 
ni = E 4aha = J2 blfiffbia = £ bib* = 1 - PJi (108) 



and to verify that : 



( ) C if3 



-I + 



E[4 



c ■ 

ia i ia 



blJiflkp + ffbiiblfi = b\M + flfi <W (109) 

= l + /7/< (no) 



effectively corresponding to the quantization of Eqs. (37) and (38). 

We can now examine, as announced above, the structure of the Fock space at site 
i Apart from the three states \fi) = fj |0) and \b ia ) = b\ a |0) , we have an infinite 
number of states of the form either || }■) = fJ(A]) n |0) or || b ia ) = b\ a (A\) n |0) , 
where (A\) n are products of b\p and fj . The constraint of Eq. (87), counting 
only the last created particle, is satisfied for all these states. The particle content of 
these states can be found using a basis of new commuting operators; for example one has 
(blMpbipbia) bUlp 1°) = b\Jlf3 1°) • Nevertheless, the states | ) and || ) share the 
same properties concerning the electron operators : 

4a \fi) = \bia) 4a IM = 4a 6 ia \fi) = 4a 6 ia \ h ip) = $a/3 |&0> 

4a II /*) = II bia) 4a II b ip) = cj Q C iQ || fi) = cJ Q C lQ || b if3 ) = 5 Q/3 || b lf3 ) 
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Thus, there is always either no electron or one electron at site i , as it is also expressed by 
the relation c| c |^ = b\ a fi b^ fi = using Eq. (102) . In spite of the infinite dimension 
of the Fock space at site i , the non double electron-occupancy at site i is well satisfied, 
and furthermore it is equivalent for physical purposes to use only the sector of the three 
states \fi) and \bi a ) . 

Let us finally insist on the fact that we have now a systematic direct algebraic proce- 
dure to find the expression in the slave-particle approach of any operator initially written 
in terms of the electron operators. For instance, since using our result fifj = 1 one has : 

= blfiffbiP = tfjip (HI) 
we obtain directly the t- J model Hamiltonian of Eq. (2) in the slave- fermion representation: 

H = t Yl (HMjho i + h.c.)+ j J2 lE&LMJ/Aa - a*^) +/*£///<- ?n 

<ij>,a <ij> a,f3 i 

(112) 

with hi given by Eq. (108) , and where we have added the \i chemical potential term, 
TV being the number of lattice sites. 

b) The slave-boson representation 

The explicit quantization in the slave-boson representation of the t-J model proceeds 
exactly in the same way, and we shall only give the following results expressed, at site i , 
by: 

kb\ = l fiJl = S aT kfl = M f = (113) 

or in other terms : 



= - b\fia Sij (115) 

We therefore directly recover in this case the expressions : 

^ = b\k ^e i ) = (s ) = 5 (116) 

which through the analysis of Section III leads to the conclusion that the expected sum 
rules for the electron spectral function and for the electron number are indeed found also 
in the slave- boson approach of the t-J model within the decoupling approximation. 
With the expression of the electron operator : 

4* = flJi (117) 

and using our results for the quantization, let us also found the exact electron number 
operator : 

«i = E 4^, = E flMh* = E = 1 - h ^ ( 118 ) 

o o a 
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= ( 1 - b\bi ) 6 tj 



fiat fj T ( $ctt + 



fir fi' 



) 6^ (114) 



6,, fl =-flb % 5. 
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f 



ia, 



to be compared with the ambiguous result using the (incorrect) naive quantization : since 
boson and fermion operators would commute in this naive quantization, c\ a c icr could 
either be written fj a fi a bibl which using the slave- particle constraint and the naive 
commutator of the boson would give ( 1 — (b\bi) 2 ) , or be written ^ a bifl a fiJb\ which 
in the same way would give — (b\bi) ( 1 + (b\bi) ) • 
Finally, since using our result hi b\ = 1 one has : 



flMfir = flfir (119) 

we obtain directly the t-J model Hamiltonian of Eq. (2) in the slave-boson representation: 



<ij>,°- 



<ij> cr,T 



(120) 

with hi given by Eq. (118) , and where we have added the fi chemical potential term, 
N being the number of lattice sites. 

V - Generalization to a slave-boson approach of the Hubbard model 

A slave-boson approach was introduced in Ref. [6] for the Hubbard model, i.e. without 
neglecting the doubly occupied sites, and it was also claimed in Ref. [5] that the sum rule 
of the electron number was still violated there. Let us show explicitly that it is not the 
case if one uses our present quantization. 

In Ref. [6] the electron operator was written in the finite- U Hubbard model as : 

cl = flb li + ab\J i{ _ a) (a = ±l) (121) 

the bosons b\ and b\ describing respectively the empty and doubly-occupied states. The 
slave particle constraint here is : 



blihi + blhi + J^flfi^^ 



(122) 



The explicit quantization of this slave-boson representation of the Hubbard model 
proceeds in the same way as shown above for the slave-particle representation of the t-J 
model, and we obtain the results at site i : 



b a i = 5 a p f ia f ri 



f =5 

IT U OT 



or in other terms : 



bah bpj 



b ■ f 



( <W - b^bai ) 5. 



bai /, 



fia i f 



t 







f b j 







(123) 



( S aT + flh, ) Sn (124) 



f b ] 



*3 



(125) 
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where a or (3 = 1, 2 and a or r = ±1 . 

The electron field being gauge invariant, one verifies as above that at the classical 
level one can use the initial brackets of the b 's and / 's to compute the Dirac brackets of 
the electron field : 

* { C L^ C ir} =i { C t> C ir}^= (126) 

while at the quantum level our results give : 

4^ = /X/<T + &£Aikr ( 127 ) 

kAo = fl_ r) f i{ _ a) + b\ i b li 8 aT (128) 

which effectively leads to : 

= S„ (129) 
One has from Eq. (127) the expression of the exact electron number operator : 

E4a, = i- (SJAi-^Ai) (iso) 

Let us define A as the average number of empty sites, d the average number of doubly 
occupied sites and 6 the hole doping concentration : 

A=(b\ i b li ) d=(b\ l b 2l ) 5 = A-d (131) 

The total electron number per site is then correctly obtained from Eq. (130) as : 

n e = (1 - 5) (132) 

Let us again emphasize that the expression of the Hamiltonian in the slave-particle 
approach is obtained, using our results, through a direct algebraic procedure. Eqs. (127) 
and (123) lead to : 

n,+n,_ = (f}J i+ + bUk,) (/7_/ 4 _ + SJAi) = bi z b 2t (133) 

and we obtain directly the Hubbard Hamiltonian of Eq. (1) with nearest- neighbor hopping 
in the slave-boson representation of Eq. (121) : 

<ij>,cr <ij> 

+ ^E^A l + ^E( st iA l -^AJ- (134) 

i i 

where we have added the fi chemical potential term, N being the number of lattice sites. 



at 
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The decoupling approximation for the Matsubara electron Green's function in imagi- 
nary time is here expressed by : 



E a (r, t) = - F a (r, r) B^-r, -r) + F_ CT (-r, -r) B 2 (r, r) 



(135) 



On one hand, one obtains the sum rule for the electron spectral function within the 
decoupling approximation : 



(136) 



/-too 

= wJ2 nb i (?) ( J2 [f(9+k)a, fl+k)*\ ) + w J2 n f^ + fc ) 

+ jrJ2 n f°(ri \ [^(i+k)A( q +k)]_) + w^2 n b2((i + k) (J2 [fi°ifi*\ + ) 

q,a ^ q \ a / 

which using our quantization gives : 

/+00 
^A ea (k,u) (137) 
-00 

= A (3 - A - d) + (1 - A - d) (1 - d) + (1 - A - d) (1 - A) + d (3 - A - d) = 2 

as it should be from Eq. (129) . 

On the other hand, Eq. (135) gives the expression for the number of electrons in state 
k within the decoupling approximation : 



n. 



' qa 



= ~ W [(f(q+k)af(q+k)<r) Q^l q) + (&2 (q+k) ^2 (q+k) ) (fqafq 

q,a 

= £l>/a(? + fc) ( + ^X)w /<T (g + fe)n 6l ( 9 ) (138) 

+ + hen Pq* /) ~ j^^2 n b 2 {q + k ) n fo{q) 



q,a 



q,a 



Using our results for the quantization of the b\ boson and the / fermion, one obtains : 
(k) = (l-A) 2 -d 2 + 2d + jjY,n f M + k)n hl (q) j? Yl n ^ (« + k ) n f°^ ( 139 ) 



n 



q,a 



q,u 



from which we indeed get the expected sum rule of the electron number : 

jf 5> e (fc) = (1 - A + d) = (1 - 8) (140) 
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This last result also directly follows from Eq. (135) : 

«e = iE( & A) (E(&/*)) + *E(E( 

= l(l-A-d) + 2d=(l-<y) (141) 
using the results of our quantization, instead of : 

(1 + A) (1 - A - d) + (2 - (1 - A - d)) d = (1 - A 2 + d 2 ) 
using the (incorrect) naive quantization of Refs. [5] and [6] . 

VI - Conclusions 

In this paper, we have first presented at the classical level the consistent Hamiltonian 
formulation of models having a slave-particle constraint for their fields. Due to this con- 
straint, the naive canonical relations are replaced by modified canonical relations which 
are compatible with the constraint. This is achieve through the use of Dirac brackets, after 
fixing the gauge generated by the slave-particle constraint. 

We have then shown at the quantum level, for the slave-fermion and the slave-boson 
representations of the t-J model and for a slave- boson representation of the Hubbard 
model, that a consistent quantization of these modified canonical relations changes the 
naive sum rules for the slave particles. These naive sum rules used in Ref. [5] were there 
shown to lead to difficulties in these slave-particle approaches for the t-J and Hubbard 
models, coming from the fact that the sum rule of the electron number was violated within 
a decoupling approximation. On the contrary, we find that, using our quantization and 
modified sum rules for the slave particles, the sum rule of the electron number is in fact 
well obeyed. On the other hand, we obtain a systematic direct algebraic procedure to find 
the exact expression in a slave-particle approach of any operator, e.g. the Hamiltonian, 
initially written in terms of the electron operators. 

We thus show that one has to be careful about the canonical relations when using 
a direct quantum operator approach for slave-particle theories. Of course we are not 
concerned for these theories neither with the functional integral approach nor with the 
Abrikosov [14] method used in some slave-particle theories [15]. 
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Appendix 

In this Appendix, we show how the brackets of Eq. (7) and the first class constraints 
of Eq. (8) are obtained and we present the expression of the corresponding Hamiltonian. 

Let us consider the classical Lagrangian (written in real time) for n bosons b a and 
m fermions f a (Grassmann variables) on a lattice : 

L = <J(A + 1) JXft&ia +*§(A - 1) 5>&U« 

i,a i,a 

+ A(&+, b, /t, /) + £ A, ( blfiia + E fUi* ~ 1 ) 

i a a 

the Lagrangian of Eq. (5) being obtained for A = A' = 1 . The canonical momenta are 
then : 

b i(x =i\{\ + \)b\ a b\ a = i\{\-l)b ia 

J ia = -q(X' + l)fl Jl = - 1) fu, n, = (A2) 
The canonical graded (see Eq. (9)) Poisson brackets, the non zero ones being : 

{fiaJ jT } = {flj] T } =-8i S 8*r {A*,n,} = (A3) 

are however not compatible with the expressions of the canonical momenta. Following the 
general procedure of Dirac [7] , one has to consider here as Hamiltonian : 

H 1= - X(&t , b, /t, /) - £ A, ( £ + £ - 1 ) 

i a u 

+ e (BiaUia + 4*4*) + E + v l F l) + E ( A4 ) 

where the u , it^, u , i^, mj are unknown (at this stage independent of the fields and of 
the canonical momenta) coefficients and where : 

B lct = l ia - + 1) b\ a Bj a = 11 - - 1) b ia n< 

F ia = J ia + + 1) f\ a Ft = Jl - - 1) f ia (A5) 

are primary constraints which are weakly zero, meaning that one has to set the constraints 
only after computing all the brackets. In order to have a consistent system, we require the 
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time derivatives of these primary constraints to be weakly zero ( ~ ), which gives : 

dX ( i dX 

{B lC(1 H 1 }=— + Aib\ a - i u\ a « [Bl, H 1 j = —^ + A t b ia + ««»~0 (A6) 

BX ( i BX 

{F ia ,H 1 } = --- Ajl + ivl^O {pl ri H 1 } = — r +A i f ia - iv ia *0(A7) 

B}ia >- 
a cr 

Eqs. (A6) and (A7) determine the u , it" 1 ", v , and Eq. (A8) gives the slave-particle 
constraint which appears as a secondary constraint. Repeating the process, we require the 
time derivative of this secondary constraint to be weakly zero, which gives : 

H 1 } = Y J (bluia + u\ a b ia ) +£(- f} a v ia + vlfur) « (A9) 

a <r 

Using the expressions of the u , -u^, v , i;t given by Eqs. (A6) and (A7), the Eq. (A9) 
reads : 



ct ice 



which is identically satisfied if X is a function of b^b and f^f . 

One recall that a quantity is called first class if its bracket with each of the (primary 
and secondary) constraints is weakly zero, and second class if at least one of these brackets 
is not weakly zero. Hi is then first class, while our other constraints are a priori second 
class since : 

|s iQ , B]^ = - i Sij S aP F] T } = - i Sij 8 aT (All) 

{B w = - bljij {Bl, *,-} = - 6 4o (5 y (A12) 

{F^ $ .} = + fl5 tJ {Fl, $, } = - /^y (A13) 



However, the following linear combination of these second class constraints : 

vjr . = * 4 + (- + + i £ (F^ + (A14) 

a a 

= i £ (- &«.&<« + + • £ (7-/- - - 1 (^15) 

(which is weakly equal to <E>i ) can be verified to be first class. 

Thus, Ilj and are first class; and B ia , Bj a , F ia and f\ g (which are such that no 
linear combination of them is first class) are second class. Defining (pi = B ia , ip 2 = Bj a 
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, (p% = F ia , (f4 = , the matrix C a b = {^a^fe} is non singular. Systematic use of the 
standard Dirac bracket [7] of two quantities A and B : 

{A, B} D = {A, B}-J2 {A Va} (C-% B} (A16) 

a,b 

then allows to set all these second class constraints strongly to zero because the Dirac 
bracket of anything with a second class constraint vanishes. 
One then finds the following results : 

{ b l> h )p] D = U X - 1) % <W = ~U X + 1) <W 

{&L^} S = ^i(A 2 - 1) <W {& ia , 6]^} d = - % S^ap (A18) 

{lL f] T } D = - 1) % <W {7i„, = -§(A' + l) <5y <W (A19) 

{7l,7 iT } D = "*l(A /2 - l)M*r {/^/ir} D = - i M«rr (^0) 

{A,,n j } D = ^ J (A21) 

Since we can now set the second class constraints strongly to zero, i.e. use the expressions 

of bi a , b\ a , f ia , f\ a given by Eqs. (A2), a choice of independent non zero canonical 
relations is : 

{b ia , b ]p} D = ~ i 8ij&a0 {fia, fjr} D = ~ i $ij$vT {A*, Kj} D = $ij (A22) 

which are those of Eq. (7) of Section II where we have drop for convenience the subscript 
D . Let us note that the parameters A and A' of the Lagrangian (Al) no longer appear, 
and can thus, as usual, be taken as unity in the Lagrangian. 
At this stage, the Hamiltonian is : 

h 2 = - x(b\ b, f\ f) - e a, ( E b l b ^ + E fit* - 1 ) + E u ^ ( A23 ) 

i ct a i 

and one now has to deal, as shown in Section II, with the first class constraints : 

*i = E b l b ™ + E fUur -1~0 Hi « (A24) 

a a 

Let us note that we have consistently for example {bi a ,Hi} = {bi a , H 2 } D since : 

{h^Hx} = u ia (A25) 
dX 

{b ia ,H 2 } D = i— 1 - + iA i bi a (A26) 
db 1 



la 



which are indeed equal using Eq. (A6) . 
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